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Abstract 



The existence of global solutions is established for compressible 
Navier-Stokes equations by taking into account the radiative and re- 
active processes, when the heat conductivity k + Q'^) < «^ < 
^2(1 + > 0), where 9 is temperature. This improves the pre- 
vious results by enlarging the scope of q including the constant heat 
conductivity. 

1 Introduction 

In this paper, we study a free boundary problem of viscous, compressible, 
radiative and reactive gases which are driven by self-gravitation. The system 
considered can be described by the conservation of mass, the conservation of 
momentum, the balance of energy and the react ion- diffusion equations 



P{et + UCy) = in9y)y + ( -]9 + fiUy)Uy + XpC^Z"^, 

^ p{zt + uzy) = {dpzy)y " p^z'", a{t) < ?/ < t > 0, 
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where p = p{y,t), u = u{y,t), z = z{y,t) denote the density, velocity, mass 
fraction of reactant of flows, respectively. In addition, the positive constants 
p, d and A are the coefficients of the viscosity, species diffusion and differ- 
ence in heat between the reactant and the product. a{t) and b{t) are free 
boundaries defined by 

a'{t) = u{a{t),t), h'it) = u{b{t),t). 

f = fiy^t) is the external force determined by / = —Uy where U solves the 
following boundary value problem: 

Uyy = Gp, a{t) <y < b{t), 
U{d,t) = 0, d = a{t), hit), 

where G is the Newtonian gravitational constant. 

Moreover, the reactant rate function cf) = (f){p, 6) is given by the Arrhenius- 
type law 

^{p,e) = Kp^-'e^exp(^-j^, (1.2) 

where m > 1, /3 is a non- negative constant, while K and A are positive 
constants. 

Finally, the pressure p and the internal energy e can be formulated as 
follows: 

p = Rpe + -e\ e = CJ + a—. (1.3) 
3 p 

where the positive C^, is the specific heat at constant volume, a > is the 
Stefan-Boltzmann constant and R is the perfect gas one. From the exper- 
iments results for gases at very high temperature, refer to ^2], the heat 
conductivity is 

Ki{l + 9'^) <n{p,9) <K2{l + 9''), q>0, (1.4) 

where Ki and 1^2 are positive constant. 

The system (II. ip is to be completed with initial-boundary conditions. 
More precisely, we suppose here that 

{p,u,e,z)\t=o = {po{y),uo{y),0o{y),zo{y)), (1.5) 
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and 

{— p + fiUy,k9y,dpZy){d,t)) = {—pe, 0,0), d = a(t),b(t). (1.6) 

The same initial-boundary problem is proposed by Umehara and Tani in 
[22l[23]. We aim at extending their results by enlarging the scope of q and (3. 

The system (11. ip has been extensively studied from the wellposedness of 
global solutions to asymptotic behavior of solutions under various initial- 
boundary conditions since the radiation is an important factor in astro- 
physics. Among them, Documet [12] considered the combined radiative and 
reactive case only for R = and then established the existence of a global- 
in-time solution for the Dirichlet boundary problem with one order kinetics 
in [13] when g > 4. In another paper [15], Ducomet and Zlotnik treated the 
case of higher order kinetics for fairly general kinetics law when q > 2. In 
addition, Ducomet and Zlotnik [TB] investigated the global-in-time bounds 
and exponential stabilization for solutions in L'^ and norms by construct- 
ing new global Lyapunov functions when q > 1. Moreover, Umehara and 
Tani [22j showed the global existence of solutions for a self-gravitating, vis- 
cous, radiative and reactive gas in the case of the external pressure Pe > 
on the free boundaries when 4 < g < 16 and < /3 < y and extended that 
result to g > 3 and < /3 < g + 9 in [23j. Qin-Hu-Wang [20] filled the gap 
with 2 < g < 3 and < /3 < 2q + 6. There are also various excellent works 
on the multi-dimensional radiation hydrodynamic equations in this direction, 
see [inmilEI] and references cited therein. 

However, those results are far from satisfactory. For example, the con- 
stant heat conductivity, which is reasonable in physics, is excluded due to the 
restriction of mathematical techniques. Our aim, in this paper, is to show 
global existence of a radiative and reactive gas for a free boundary prob- 
lem including constant heat-conductivity (g = in (II. 4p ). We remark that 
similar problems for perfect or real fluid flows have been extensively studied 
in [T|-[8| ITn - [T9| [2l] and references cited therein by many researchers. However, 
those methods adopted in the literature is not valid for the case of radiative 
flows, which arise extra difficulties in mathematical analysis. More precisely, 
besides the possible concentration of mass and heat, the radiation equation 
and (p, 9) dependence of heat conductivity make it comphcated. 

More recently, Zhang-Zhang [25] established global existence result for the 
system (II. ip to Dirichlet boundary problem with respect to velocity when 
g = and (3 G [0,8] . However, it leaves a gap for the free boundary 
boundary problem suggested in [221 123]. We are intend to make the gap 
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in this paper. The novelty of this paper is to provide a new insight into 
solving radiative and reactive flows. Following the framework of [22] in the 
context, the main obstacle of a priori estimate is for Vx € L°°(0,T; L^[0, 1]), 
see Lemma [3T8l which in turn is coupled with the integrability of temperature, 
and the "regularity" of heat conductivity. To overcome the above difficulties, 
we introduce the auxiliary function G{v, 6), refer to Lemma [3.61 as in P^[T^ 
and succeed in getting ^ G L'"(0, T; L°°[0, 1]), < r < g + 5, g > 0. Then, 
using the standard estimation of parabolic equation, we further get 6 G 
^4n-4j-g^ l] wheu q > and < /3 < q + 9. With those essential a priori 
estimates, we can deduce a priori estimates of higher derivatives of {v, 9, u, z). 
This methods not only adapt to the Dirichlet problem in [2S] but also may 
be shed light on the global existence of various radiative models, which is 
our further study. 

The outline of the paper is as follows. In section 2, we introduce the 
Lagrangian mass coordinates to transform the free boundary problem (11. ip . 
flL5l) - fll.6l) into the fixed boundary problem fl2.ip - fl2.3p . which are equivalent 
each other. In section 3, we deduce some a priori estimates for state param- 
eters (p, M, 6*, z), for which it is necessary to guarantee the global existence of 
solutions by extending the local solutions. 



2 Main Result 



In order to study a free boundary problem conveniently, we introduce the 
Lagrangian masss coordinates transformation. Let 

x = r p{i,t)di, t=t. 

Ja{t) 

Then, the free boundaries a{t) and h{t) become a; = and x = J^^^^^ p(^, t)d^ = 
M, which denotes the total mass of fiuid, respectively. By the conservation 
law of mass, we may assume M = J^j^j^^ p{^,t)d^ = J^j^^^ p{^,0)d^ = 1, i.e., 
< X < 1. So the system (11. ip becomes in Lagrangian mass coordinates 
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vt = Ux, (2.1a) 
«t=f-p+— ) -g(x-I), (2.1b) 



et=(-e^ +(-p+^)«,. + A#-, (2.1c) 
= a^-) - "f"^""^ < X < 1, t > 0. (2.1d) 



V 

The corresponding boundary conditions are 



= (-Pe,0,0), (2.2) 

x=0,l 



and the initial conditions are 

{v,u,9,z)\t=o = ivo{x),Uo{x),9o{x),zo{x)), xG[0,1]. (2.3) 

where f = 1/p is the specific volume and u—J^ udx is instead of u. We remark 
that the external force / is more realistic than that in [H] by Ducomet. 
The main result is stated as follows. 

Theorem 2.1. Let ffL2D-ffLD hold, and suppose that {q, f3) E [0, +oo) x 
[0,g + 9),Pe > 0. Moreover, the heat- conductivity k in (11.41) satisfies 

\Kp{p,d)\<K2{i + ey, \ne{p,9)\<K2{i + 9y-\ 

Assume that 

(t;o(x),no(x),^o(x),zo(x)) G C^+"([0, 1]) x (^^+"([0, 1]))^, 

and 

vo{x), e^{x) > 0, < zo{x) < 1, xe [0, 1], 

where a G (0, 1). Then there exists a unique solution (f , m, 9, z) of the initial 
boundary problem fl2.1l) - fl2.3p for any fixed T > 0, such that 

{v, vx, Vt) G (C"'"/^(QT))^ («, z) G (C^+"'^+"/2(QT))^ 
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and 

v{x,t), 6{x,t)>0, 0<z{x,t)<l, X eQ^, 

where Qt = (0,1) x (0,T). 

Here, as usual, C™"'""[0, 1] denotes the Holder space on [0,1] with expo- 
nents me N+,ae (0, 1), C"'"/2(Q^) the Hdlder space on Qt = [0, 1] x [0, T] 
with exponent a in x and a/2 in t with a G (0, 1). Similarly, if a function 
f{x,t) e C2+a,iW2 ^^th a G (0, 1), then it means K^,Mt) e C"'"/2(Qt)- 

Remark 2.1. The heat conductivity in (ll.4p covers the case of 

by K = Ki + K2J. 

Remark 2.2. The external force p^ > is not essential. The result are also 
satisfactory for pe G R provided that we get v G L^[0, l],see Appendix below. 

The existence and unique proof of Theorem 12 .11 is motivated by the proce- 
dure devised in [9] through the apphcation of the Leray-Schauder fixed point 
theorem, see also [9l[T8]. Therefore, a priori estimates of the local solution 
{p,u,9,z) and its derivatives are necessary. 



3 Proof of Theorem 2.1 



In this section, we establish some a priori estimates for the state parameters 
(p, u, 9, z) so as to get the global existence of solutions to the problem (12. ip - 
(12. 3p . In the sequel, the generic positive constant C{C{T)) may be different 
from line to line. In addition, we denote := sup^^^^gQ^ |m(x, t)| and 

/3+ = max{/3,0}. 

The following a priori estimates are basic and useful. Please refer to [22] 
for detailed proofs. 

Lemma 3.1. Under the hypotheses of Theorem \2.1l it satisfies that 

he + Xz-{ ^ —\dx + pe vdx<C, (3.1) 

V 2 2 / Jo 

U{t) + / V{s)ds < C{T), 

and 

< z{x,t) < 1, 







^'^dx+ r I'^zld-'" ' ^* f ^~n^+l^^^.._ f^4i^) 



-dx + / —z^dxds + (pz'^^^dxds = / dx, 

iQ ^ Jo Jo Jo Jo Jo 2 
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where 

U{t) 

and 



[ [C^{9 - 1 -\og9) + R{v -l-\ogv)]dx, 
Jo 



forO<t<T. 

In the sequence, it gives the upper bound of density. 
Lemma 3.2. There holds that 

C,{T)<v{x,t)^pix,t) < (3.2) 
for < X < 1 and <t <T. 

Proof. Integrating fl2.1bp over [0,a;] x [0,t], we obtain 

1 r. N 1 1 1 /"* , Gx{l-x)t pe 
- {u-Uo)d^ = logv -logvo / pdrH \ 1, 

that is 

1 /•* , Gx(l-X)t 1 r\ ^ Pe , , 

logt; = logt;o + - / pdr \- ^ + - / {u-uo)dC--t. 3.3 

Thus, we find by f l3.1jl that 

I r , , ,^ Gx(l - X)t Pe 

log > log Wo + - / [u- Uo)d^ 1 

>logfo-c/ u^dx - -( ^ + pe]t - C 



JO / 
> \ogvo - Ct - C, 

which leads to ([32]). □ 
The following auxiliary result stems from the idea of |22] . 
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Lemma 3.3. One has 



max6''"(x, < C(T), 



where < r < g + 4 and g > 0. 

Proof. It is known from the mean value theorem 



9dx, 



where x{t) G [0, 1] for every t G [0, T]. Furthermore, we arrive at 



^ Jxit) 



< C7 1 + 



v^Q 



< C 



< C 



1+ ( l^—dx ] V{t)i 



1+ / {v + ve*)dx] 



(3.4) 



due to0<r<g + 4 and Young inequahty. Thus it is complete by integrating 
with respect to time. □ 



Similarly, we obtain the lower bound of density, which is very important 
to obtain the existence of global solutions. 

Lemma 3.4. If q>0, then it holds that 

v{x,t) < C{T) ^ p{x,t) > 1/C(T), 
for {x,t) G (0,1) X (0,T). 
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Proof. It follows from fl3.3p that 

1 /■* , Gx(l-x)t 1 f\ , pe 
logv = logvo + - / pdr h-/ (u-Mo)d4 1 

< logt;o + C /" max(^ + ^^)cis 
Jo [0.1] 

+ C /" u^rfx + C /" uldx + C{T) 
Jo Jo 

<\ogvo + C{T), 

by Lemmas 13.21 and 13.31 Thus it is complete. □ 
Lemma 3.5. Similarly, we have when q > 

[ [ uldxds <C{T). (3.5) 
Jo Jo 

Proof. The equality ( ]2.1b[l implies by integrating from (0, 1) x (0,t) that 



^ [u^ Gx{l — x)v 



dx + n / —dxds 
'o Jo V 



^ [u^ Gx{l — x)v 



{x,Q)dx + Pe / VQ{x)dx 



V 2 2 y Jo 

+ / / pUxdxds — Pe / vdxds 
Jo Jo Jo Jo 



<C{T) + e I I uldxds + / / p'^dxds, 
Jo Jo Jo Jo 

which implies for fixed small e that 

^ ^ uldxds < C{T) + ^ max^^ (^j^ O^dx^ ds 
< C{T). 

The proof is complete. □ 
Lemma 3.6. For < £ < 1 and q > 0, there holds 



^0 



Ql+e 



dxds < C. (3.6) 



Proof. Let G{v,9) = ^-'e^{v, ^d^. Then 

pO 2 

Integrating it over [0, 1] x [0,t], we arrive at 



Jo "»'+■ Jo 7o Vi'f !}• 

Gdx — / Goc^a; + {1 — e) / / ^ "^Ps^d^ ■ Uxdxds 



^0 ^0 ^0 

t rl i-t nl 



Jo io Jo Jo 

< c. 

□ 

With the proceeding a priori estimates at hand, we can firstly improve 
the integrabihty of temperature and deduce that 

6 e L"(0, T; L°°[0, 1]), < r < g + 5, g > 0, (3.7) 

when the function V{t) = jrf^dx is instead of V{t) in the proof of Lemma 
I3.3[ Furthermore, we can make a priori estimates of higher integrabihty of 

e. 

Lemma 3.7. If0</3<q + 9 and q > 0, then it has 

[ 9^''-^dx+ [ max^«+"c/s+ / / {\ux\'' + k9'^''-^91) dxds < Cn, 
Jo Jo Jo Jo 

(3.8) 

for sufficiently large n > 2. 

Proof Multiplying e^""-^ on Olcj) . we get 

1 ft pi n4n-9n2 

9'^-^dx+ / / -^dxds 



^0 



<c + c [ [ {e^'^'^ul + e^""-^ \ux\ + 6^+^''-'') dxds 

Jo Jo 

<C + Cn [ [ \uTdxds + Cn! max(9^ + 9^^'^'^+) [ O^^'-^dxds. 
Jo Jo Jo 10,1] 7o 



(3.9) 
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On the other hand, we can rewrite (12.1bp with the initial boundary data 

II I Gx(l—x) I 
Wt = ^Wcc - p H 5^ '- + Pe, 

w\t=o = woix), 
W\x=0,l = 0, 

where w = ud^. The standard estimates of solutions to linear parabolic 
problem give rise to 

ft ^1 



"f 

^0 



Ux\^dxds 



'0 ^0 

<c(i- 



Wxx.rdxds 

p^dxds 



< C 1 



Jo 

t rl 



e^'^dxds 



< C 1 



Jo 

t 



which together with (13.91) leads to 



t /.I 



Jo 



max^^ / e^^'-^dxds I . 

[0,1] 



K9"etdxds 



<C + C [ max(9^ + 9^^-"^+) [ 9^''-^dxds, 
Jo [0.1] Jo 

which implies 9 G L^"~'^[0, 1] by Gronwall's inequality which in turn, it gives 
Ux G ((0, 1) X (0,t)) according to (IHl) and in the same manner, it implies 
that 9 e L«+^'^-4(0, T; L°°[0, 1]). Thus, it finishes the proof. □ 

Lemma 3.8. If0<f3<q + 9 and q > 0, we arrive at 



v^dx 



JO 



9v,.dxds < C, 



Proof. Multiplying ( I2.1bl) by m — and integrating it over (0, 1), it yields 



1 d 



~~r ^x dx + 







2dt 

R9uu 

y2 



v^dx 



■dx — 



(^u — —Vx^ dx. 

(3.10) 
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On one hand, we deduce that 



^ ROuv, 



dx 



<s Ov^dx + CsTonaxO ■ / u dx 

[0,1] 







(3.11) 







<e I Ovldx + Climax 9 (x,t). 

~ Jo m 



On the other hand, we get 



^ + fe']e^ + G{x 



< c 



< 



1+ [ %dx+ [ e^eldx+ [ {1 + e^) (u - ^vX dx 

Jo ^ Jo Jo \ V / 

C 1+ -^dx+ / K9^''-^9ldx+ / {l + e^){u^ + vl)dx 
Jo "^^ Jo Jo 



< C 



1 + / '^dx + / KO^'^-^eldx + ( 1 + max^' 1 ■ ( 1 + / vidx 



[0,1] 



which together with f l3.10p and (13. lip for any fixed small e leads to 



/ vidx +11 
Jo Jo Jo 



Ovtdxds 

Jo 



\Jo 



<C + C I ( / -^dx+ I Ke"eidx + max (x, s) ] ds 



[0,1] 



+ C I ( 1 + max02 

[0,1] 



<C + C I ( l + maxO^ 

[0,1] 



v^dx I ds 



v^dx I ds 



Thus, the lemma is completed by Gronwall's inequality. 



□ 
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We introduce the following new variables 

X := [ [ {l + e'')e'^dxds, (3.12) 



^0 



Y := max [ {1 + 9^'')9ldx, (3.13) 

0<t<T Jq 

Z := max / ul^dx. (3.14) 

0<i<T /n 



Firstly, with the aid of interpolation and embedding theorem, we deduce 
that 



\4'^<c{i + zl) 



by simple calculation. 

Lemma 3.9. Under the assumptions of Theorem \2.1\ we have 

max eix,t) < C + CF2u+2n-i), (3.15) 

[0,l]x[0,i] 

for {x,t) G (0,1) X (0,T). 

Proof. By the embedding theorem, we deduce that 

meix9''+^"-\x,t)<C [ O'^+^^'-^dx + Cn f (1 + e)^+'""'|^.Ma;, 

[0,1] Jo io 

which implies that 
max0'^+2"-i(x,t) 

[0,1] 

< CmaX^''+2n-2 I" Qd^j^Cnl (1 + ^)'^+" I^J rfx 
[0,1] 



< emax^'?+2"-i + C„ (^j^ (1 + 9fWldx^ ' (1 + ^)^"-^dx^ ' + C,. 

So we deduce f l3.15p . □ 
Lemma 3.10. One has 

X + Y <C (l + (3.16) 
for < (3 < q + 9 and q>0. 
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Proof. We introduce the function as in 

Jo V 

The simple calculation leads to 

Kt = -Ot + K,u,, (3.17) 

V 

K^t=(-6x) + K^vVxUx + K^u^^+ (-) vJt, (3.18) 

\V Jt \V/v 



\K,l\K,,\<C{l + ef^\ (3.19) 
Multiplying f l2.1cp by Kt and integrating it over (0, 1) x (0,/!:), we get 



t pi n2 pt /•! 





'0 ^0 ^ Jo Jo V 

K^Uxdxdr — { OvaUx m?. ) —OtdxdT 

Jn Jn V V V 



'0 Jo Jo Jo 



- / / ( - ) K^Uxdxdr - / -O^K^vVxUxdxdT (3.20) 
Jo Jo ^ V ^ Jo Jo V 

- / -OxK^Uxxdxdr - / -6*^. ( - ) VxOtdxdr 
Jo Jo V Jo Jo V \v/v 



+ \ [ [ (l)z^KtdxdT. 
Jo Jo 

In the sequel, we will estimate all terms in f l3.20p . 
Firstly, we have 

t rl Ci2 rt rl 

'-dxdr>c I {i + e^){i + e'i)eldxdT>cx, (3.21) 




and 



JO ^ JO JO 



I [ -ex (-ex) dxdr 

Jo Jo V \v )t 

>c I (1 + e'ifeldx -c>CY -c, 
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by the definition of flXT^ and (TO]) . 
Secondly, we find that 



t /•! 




egOtK^Uxdxdr 



^0 



t rl 



< c 




^0 




u^dxdr 



JQ 



<eX + C'eF2(, + 2n-l) 

< e{X + Y) + C, 
Similarly, we have 

j j [dpeUx - ^ul^ '^dtdxdr 

<c{T) [ [ [{1 + ey+^\u,et\ + {1 + ey\et\ul] dxdr 

Jo Jo 



<eX + Cs\{l + 6'^+^) 



|(0) 



t 




u^dxdr 



JO 



< e{X + Y) 



91(0) 



t fl 




u,.dxdT 



JO 



and 



J J {dpeUx - ^ul^ KyUxdxdr 

<c [ [ [(1 + ey+^ul + (1 + 9y+^ \uxf] dxdr 

Jo Jo 



h5M(0) 



t p1 




u^dxdr 



<eY + C^ 



g+l|(0) 



JO 
t rl 




\uJ dxdr + C 



JO 
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Likewise, we get 



^0 



—OxKyyVxUxdxdr 

V 



< 



C f j {l + ef''^^%VxUx\dxdT 
Jo Jo 

i7'i"^)'(i M<^+ 



< C \Ux\ 

< CZi. 

In the same manner, we find that 
-t /.I 




K 



JO 



V 



^ '^'XJj r^rj^ dX 



At last, one has 

ft rl 




—dj. (—] VxOfdxdr 

lo Jo ^vJv 



< 



c f f {i + ef^iexVxetidxdr 

Jo Jo 

2 / --i (1 + efivy 



<eX + a I max( — 



[0.1] V V 



dx I dr 



V 



vidx I dr 



max I dr. 

[0,1] V V 



(3.23) 



<£X + CJ(l + e)^|(°) / max( - 

Jo [0.1] V 

Jo 

The next goal is to show the estimation of the last term on the right hand 
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side of the above inequality. From the embedding theorem, it yields 



[0,1] \ V J 



V \ V 



dx dr. 



9+2 

< (7^2(9+2,1-1) 



g+2 

_|_ 4{9 + 2n-l) 




JO 



—^dxdr 
v9^ 



V 



9 + 2 9 + 2 

< (7F2(9 + 2n-l) _|_ 4(9+2n-l) 



(3.24) 

In the following, we derive the estimations of all terms of the right hand 
side in ^2^. 

Firstly, from f l3.12p and flS.lSp . we get 

f [ {i + ey+^eleldxdT 

Jo Jo 

<c|(i + 0)Y°^/ / {i + efeldxdr 

Jo Jo 

<cUi + eff^x 



< C (^X + r2(9+2n-l)X 

Similarly, recalling (13. 5p and (I3.15p . we deduce that 

r f {i + ey+^eYeuldxdr 

Jo Jo 



Jo 

<C [ [ {1 + 9y^^%ldxdT 
Jo Jo 



/ 9 + 10 \ 

< C ( 1 + F2(9 + 2n-l) j 
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and 



/ / {1 + ey-'^ldxdT 

Jo Jo 



-2|(0) / / „4 



<c ui + ey+Y / / utdxdr 




JO 



/ q + 2 \ 

< C M +r2U+2n-l) j . 



At the same time, we obtain 

/* [\l + 9y^''c^''z^"'dxdT 
Jo Jo 

9+2+2/3 

< (7y2(,+2n-i) + c, 



Thus 



y2(g+2n-i) / Hiax -6'^. | c/r 
/o [0.1] 



q+l 3q+2 / 8 1 1 9+10 g+2+2;3 \ 

< -J+an-l _(_ (7y4(9 + 2n-l) I _|_ y 4(<7 + 2n- 1) 2 -|- ^2 + y4(9 + 2n-l) -)- y4(g+2n-l) j 

< e{X + Y)+C, 

which together with ( 13.23^ leads to 




JO 



( — 1 VxOtdxdr 

V 



< e{X + Y) + C. (3.25) 



Finally, the term A (f)z"^ Ktdxdr gives 



A /" /" (j)z^KtdxdT 
Jo Jo 

<c [ [ (t)z"'k\et\dxdT + c f [ e^+^\ux\(i)z"'dxdT 

Jo Jo Jo Jo 



<eX+\{l + ey+ff^ [ I <Pz^'^dxdT + C\d'^+f'+'\'^'^ [ I \ux\dxdT 

Jo Jo Jo Jo 

q + fi + l 

<eX + CY2(q+2n-i) 
<e{X + Y) + C. 



(3.26) 
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Collecting (K2U\f -( Kmf . we deduce flXTB]) . □ 
Lemma 3.11. We have 

f [< + el + ul, + ul) dx + /* [el + ul,) dxdr < C{T), 
Jo Jo Jo 

and 

l«^|(o) + l^|(o) + l^|(o)<c(r). 

Proof. We differentiate fl2.1bl) by t and then multiply the equality by Uf to 
obtain by integrating it over [0, 1], 



— [ —dx + n f —dx 
dt ./n 2 ./n V 



(ptUxt + ^ulu^t^ dx. 



Integrating it over [0, t] and using the Young inequality, we get 



ufdx + / u1.^dxds 
Jo Jo 



< 



< 



C + C [ [ {p^t+ul)dxdT 
Jo Jo 

C + C [ [ {l + 9y^dxdT + C\e^f^ [ [ uldxdr 
Jo Jo Jo Jo 



"t pi 
+ C I / utdxdr. 




JO 



and hence 



u^dx + / u^fdxds 
Jo Jo 

<C (l + X + Y^+^X + Y^- 



1 „ g+2n + 2 
< C 1 + ^2""^+^ 
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On the other hand, the momentum equation (12.1b|) imphes 

/ ul^dx < C + C {uf + pI + ulvl)dx 
Jo Jo 

(^1 + u'.dx + + Oldx + (\0'\^'^ + \ul\^'^) j^' vldx 



< c 



/ 1 ^ g + 2n+2 3 \ 

<c(i + z^ 



It follows from this that 

z<c(i + ziy 

which implies the boundedness of Z by Young inequality when g > and 
n is sufficiently large. So are X and Y by (13.161) . Furthermore, it holds for 
f\ W\ loi^l + + ul,)dx and J^iO^ + u^dxdr. so is u by the 



embedding theorem. This is complete. □ 

With the aid of the above a priori estimates, we further deduce a pri- 
ori estimates of higher derivatives of (p, u, 6, z) which is analogous to the 
reference [22], and then we omit it proofs for brevity. 



Lemma 3.12. It satisfies when < /9 < g + 9 and g > 0. 

0{x,t) > C, 

/ {^1 + ^Ix + \ / zl^dxdr < C, 

^0 Jo Jo 

I {Olx + olM^ +11 eldxdr < C, 
Jo Jo Jo 

For the Holder estimation of global solutions, we refer to [22] for detailed 
processes. Thus, the proof of Theorem 12.11 is complete. 

As mentioned in Remark 12.21 the same conclusion are also satisfied for 
the external force pe < on the free boundaries provided that we obtain a 
priori estimates of specific volume, see the Appendix. With Lemma 13.131 in 
hand, we can get the same a priori estimates of (f , m, 6, z). 
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Appendix 

Lemma 3.13. If Pe < 0, then it has 



1 

v{x,t)dx < C{T), 







forO <t<T. 

Proof. It is derived from integrating (I2.1b|) over [0, x] 



which imphes by multiplying v and integrating with respect to x and t that 

/i I vdx H /" f x(l — x)vdxds 

Jo 2 Jo io 

= yU / forfx — Pe / vdxds (3.27) 
+ / / vpdxds + [ [ v( [ Utdy^ dxds. 



Jo Jo Jo \Jo 



Obviously, it yields 



t .1 

/ vpdxds 
Jo 

f [Re + -voA dxds 
Jo \ 3 / 
-t pi 
< C I / edxds 



<C{T)- Cpe [ [ vdxds. 
Jo Jo 

On the other hand, we get from fl2.1al) by integration 

v{x,t) = vo{x) + / Ux{x,s)ds, 



^0 



(3.28) 
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which leads to 



/ v[ I Ufdy ]dxds 
Jo 
1 



Vo{x) + J Ux{x,s)ds j y^J udyjdx 

+ / u^dxds — I Vq{x){ I Uo{y)dy\dx 
Jo Jo Jo \Jo J 

voix)i / {u — uo{y))dy ] dx + / / u^dxds 



\J0 / JO JO 

1 







— / u{x,t)[ / u{x,s)ds\dx 



by integrating by parts. Using Young inequahty and initial conditions, we 
get from (13. ip 



J ^(^J '^tdy^ dxds 

1 rt /•! 



<e [ u^dx + Ce I [ u^dxds + Ce (3.29) 
Jo Jo Jo 



^1 l*t l*\ 

< —epe I vdx — CsPe / / vdxds + C£{T). 



10 Jo Jo 

Collecting we have 



(yU + epe) I vdx + Y / / ^(1 ~ x)vdxds 



JO 



<-CePe [ [ Vdxds + C,{T). 

Jo Jo 



Choosing sufficiently small e and by Gronwall inequality, one obtains the 
desired result for p^ < 0. The proof is complete. □ 
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